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NOTATION 


alt  Altitude 

6  Wing  span 

c  Mean  aerodynamic  chord 

Cp  Drag  force  coefficient  = 

CG  Centre  of  gravity 

Ci  Lift  force  coefficient  = 

Cr  Total  force  coefficient  =  \JC^  +  Cq 

Cx  Axial  force  coefficient  = 

Cy  Side  force  coefficient  = 

Ci  Normal  force  coefficient  = 

C,  Rolling  moment  coefficient  =  ssss±^^s£»oL JlSsUlL 

Clf,C,t,Clr  Rolling  moment  coefficient  w.r.t.  roll,  pitch,  and  yaw  rates 
respectively 

Cm  Pitching  moment  coefficient  =  SiUHsL  asrms^ ...»/  tasiS it 

Cm,, Cm,, Cm,  Pitching  moment  coefficient  w.r.t.  roll,  pitch,  and  yaw  rates 
respectively 

C„  Yawing  moment  coefficient  =  gfflsat  msssU^L^t sLessHil 

Cn, ,  Cn, ,  C„,  Yawing  moment  coefficient  w.r.t.  roll,  pitch,  and  yaw  rates 
respectively 

g  Gravitational  constant 

lx,  Iy,  lz  Moments  and  product  of  inertia 

Lgw  Vector  transformation  matrix  from  wind  axes  to  body  axes 

m  Aircraft  mass 

p  Roll  rate  about  body  axes 

q  Pitch  rate  about  body  axes 

R$  Spin  radius  measured  from  the  spin  axis  to  the  aircraft’s  CG 

r  Yaw  rate  about  body  axes 

5  Aircraft  wing  area 

V  Flight  path  velocity 
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INTRODUCTION 


This  report  documents  the  development  of  a  flight  dynamic  model  of  aircraft  spin¬ 
ning  which  has  been  written  using  the  Advanced  Continuous  Simulation  Language 
(ACSL).  The  flight  dynamic  model  has  the  facility  to  simulate  time  history  ma¬ 
noeuvres  from  both  rectilinear  flight  and  from  steady  spin  conditions  and  provides 
information  on  the  aerodynamic  forces  and  aircraft  motion.  The  model  is  applied 
to  a  basic  training  aircraft  design  which  is  intended  to  have  good  characteristics  for 
demonstrating  spin  entry,  steady  spins,  and  spin  recovery. 

The  aerodynamic  data  required  by  the  flight  dynamic  model  has  been  stored 
in  data  base  format  and  consists  of  wind  tunnel  rotary  balance  data  from  a  basic 
training  aircraft  and  computed  disturbance  data.  The  problems  with  producing 
such  a  data  base  to  cover  a  large  flight  regime  are  discussed  in  Section  3.  The  flight 
dynamic  trimming  methods  are  described  in  Section  4,  and  Section  5  presents  the 
development  of  angular  rate  equations  needed  to  calculate  the  disturbance  moment 
contributions. 

Two  example  applications  of  the  flight  dynamic  model  are  presented  in  Section  6 
to  illustrate  the  capabilities  of  the  model.  One  gives  an  example  of  equilibrium  spin 
prediction;  the  other  demonstrates  a  spin  entry  and  recovery  manoeuvre  for  which 
various  methods  of  presentation  of  the  output  time  history  data  have  been  supplied. 


2  FLIGHT  DYNAMIC  SPIN  MODEL 

The  six  degree  of  freedom  air-path  (SDOFAP)  aircraft  simulation  program  SD- 
OFAP  of  Reference  1  has  been  used  as  the  framework  for  the  Sight  dynamic  spin 
model  program  (FDSM).  This  program  is  written  in  the  Advanced  Continous  Simu¬ 
lation  Language  (ACSL)  (Reference  2)  and  can  be  supplemented  with  user-supplied 
FORTRAN  subroutines.  The  SDOFAP  model  provides  a  comprehensive  means  by 
which  to  represent  the  aircraft’s  characteristics  by  supplying  user-defined  subrou¬ 
tines  for  the  X,Y,Z,L,M,N  aerodynamic  and  propulsive  forces  and  moments.  For 
the  flight  dynamic  spin  model  the  aerodynamic  forces  and  moments  are  derived 
from  wind-tunnel  rotary-balance  data  in  conjunction  with  computed  disturbance 
moment  data. 

In  order  to  eliminate  the  problems  of  singularities  using  Euler  angle  rate  equa¬ 
tions,  the  simulation  program  employs  quaternions  to  represent  the  aircrafts  atti¬ 
tude.  Quaternion  time  derivatives  are  always  finite  and  continuous  and  only  require 
four  integrations  with  one  constraint  equation,  whereas  the  direction  cosines  gener¬ 
ated  from  the  Euler  angles  require  nine  integrations  and  six  constraint  equations. 

The  state  variables  for  the  flight  dynamic  model  that  are  integrated  over  the 
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time  history  are: 


a 

0 

V 

P 

q 

r 

n> 

n 

Tl 

X 

y 

alt 


Sight  path  variables 


body  rotation  rates 


attitude  quaternions 


displacements 


e 


Euler 

angles 


The  aerodynamic  force  and  moment  contributions  for  the  Wamira  spin  model  data 
comprise  the  following  rotary  balance  and  disturbance  terms: 


X 

— 

Xr 

Y 

= 

Yr 

Z 

= 

Zr 

L 

= 

Lr  + 

La 

M 

= 

Mr  + 

M0 

N 

= 

Nr  + 

No 

The  decomposition  of  the  aircraft  motion  into  steady  rotating  and  disturbance  com¬ 
ponents  is  described  in  Section  5.  The  control  inputs  used  to  manoeuvre  the  aircraft 
are: 

6. 

Sr 

6a 


3  SPIN  MODEL  DATA  BASE  DETAILS 

3.1  Methods  of  Aircraft  Simulation 

An  important  consideration  in  constructing  a  flight  dynamic  model  is  the  represen¬ 
tation  of  the  aerodynamic  forces  and  moments.  Reference  3  provides  a  summary  of 
dynamic  modelling  procedures  under  consideration  for  high  rotation  raie  non-linear 
manoeuvres.  These  are  listed  below: 
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1.  Polynomial  Model:  A  Taylor  expansion  of  the  force  functions  about  an 
equilibrium  condition  yielding  aerodynamic  forces  in  derivative  form. 

2.  Tobak  and  Schiff  Formulation:  A  comprehensive  development  of  mathe¬ 
matical  models  which  can  represent  all  types  of  non-linear  time  depen¬ 
dent  motion. 

3.  Data  Thble  Representation:  Aerodynamic  data  stored  in  table  form.  This 
enables  the  forces  to  be  represented  as  arbitrary  non-linear  functions  of 
a  number  of  independent  variables. 

4.  Rotary  Balance  Data:  Rotary  balances  enable  the  model  to  be  rotated 
at  constant  incidence  angles  to  the  velocity  vector,  and  over  a  range  of 
rotation  rates  creating  an  extensive  non-linear  steady  rotation  data  base 
of  information. 

Of  these  four  methods  the  Rotary  Balance  (RB)  method  was  implemented,  using 
an  extensive  RB  data  base  collected  on  a  jtli  scale  model  of  the  Wamira.  This 
experimental  data  is  augmented  by  computed  aerodynamic  moment  data  required 
to  describe  the  additional  forces  occurring  during  the  non-steady  parts  of  the  ma¬ 
noeuvres. 

3.2  Source  of  Rotary  Balance  Data  Base 

A  Wamira  spin-model  RB  data  base  has  been  developed  from  RB  spin-model  mea¬ 
surements  carried  out  at  the  NASA  Langley  Research  Centre  (LRC),  in  June  1983  by 
Bihrle  Applied  Research  Inc.  (BAR).  The  results  of  this  testing  have  been  reported 
in  References  4  and  5  and  subsequently  presented  in  graphical  form  in  Reference  6. 

Simulating  an  aircraft’s  dynamic  behaviour  using  an  aerodynamic  data  base 
requires  that  data  be  generated  to  cover  the  full  range  of  state  and  control  variables 
involved.  The  data  produced  at  NASA  LRC  covered  a  wide  range  of  angle  of  attack 
(a),  sideslip  angle  (/9),  and  spin  rate  (O)  conditions,  combined  with  variations  in 
elevator  (5,),  rudder  ($r),  and  aileron  (5a)  deflection.  For  a  given  control  deflection, 
the  a  range  was  0°  to  90°,  the  0  range  was  -10°  to  0°,  and  the  non-dimensional 
spin  rate  (|j£)  range  was  -0.6  to  0.6.  The  settings  of  control  deflection  tested  were 
fi.  :  -25°,0°,23°;6r  =  -25°,0°,25°;«„  =  -18°,0°,18°.  However  not  all  possible 
control  deflection  combinations  were  tested. 

3.3  Limits  on  Available  Data 

The  flight  dynamic  spin  model  requires  the  aerodynamic  forces  and  moments  X,Y,Z, 
L,M,N  to  be  calculated,  therefore  the  RB  data  needed  are  Cx,Cy,Cz  (axial, 
side, normal-force  non-dimensional  coefficients  respectively),  and  Ci,Cm,Cn  (roll, 
pitch,  yaw-moment  coefficients  respectively).  Assuming  that  all  the  control  inputs 
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are  coupled,  (ie.  a  change  in  rudder  deflection  will  affect  the  influence  of  the  ele¬ 
vator  etc.),  then  the  six  force  and  moment  coefficients  will  be  coupled  functions  of 
the  three  control  inputs  and  the  three  state  variables  as  follows: 

Cx  =  fn{S.,6r,$a,p,a,§}  (1) 

and  similarly  for  Cy ,Cz,Ci,Cm,C„.  The  generation  of  a  data  base  needed  to 
describe  the  non-linear  variation  of  the  six  aerodynamic  coefficients  with  the  three 
control  settings  and  the  three  flow-angle  states  allowing  for  full  coupling  would 
require  a  prohibitively  large  test  program.  For  the  Wamira,  measurements  were 
made  at  ten  angles-of-attack ,  eleven  rotation  rates,  and  two  sideslip  angles  for  a 
large  combination  of  control  deflections  with  zero  and  maximum  setting  angles.  This 
means  that  data  is  required  for  at  least  27  control  deflection  states  (ie.  combinations 
of  6,  :  -ve,  zero,  +ve  Sr  :  -ve,zero,+ve  S„  :  -tie,  zero,  +  tie)  to  provide  an  efficient 
interpolation  scheme  as  shown  in  Figure  1.  However  the  data  available  does  not 
cover  the  full  range  of  aileron  deflections  required  to  set  up  such  a  scheme. 

3.3.1  Decoupling  of  Aileron  Data 

Since  only  a  small  number  of  aileron  control  configurations  were  tested  in  Ref¬ 
erence  4,  the  possibility  of  decoupling  the  effect  of  the  aileron  from  both  the  el¬ 
evator  and  rudder  was  investigated.  By  examining  tests  carried  out  on  model 
configurations  with  constant  rudder  and  elevator,  but  with  variations  of  aileron 
(5tt  :  -18°,0”,18”),  the  degree  of  coupling  can  be  found.  The  tests  available  were: 


1. 

St 

=  -25°,  Sr 

= 

0”, 

St 

= 

0‘ 

2. 

6, 

=  -25”,  Sr 

0°, 

So 

= 

18' 

3. 

6. 

=  -25”,  Sr 

- 

-25”, 

So 

— 

-18' 

4. 

«. 

c* 

1 

II 

= 

-25”, 

So 

— 

0' 

5. 

s. 

=  -25”,  Sr 

= 

-25”, 

So 

— 

18' 

6. 

6. 

II 

o 

O 

4* 

: 

-25”, 

So 

— 

0' 

7. 

s. 

o 

O 

II 

= 

-25”, 

So 

= 

18' 

8. 

S, 

II 

o 

o 

4* 

— 

25”, 

So 

= 

-18' 

9. 

S, 

II 

o 

o 

-S» 

= 

25”, 

So 

= 

O' 

Calculating  the  difference  between  these  four  groups  of  test  data  gives  two  in¬ 
cremental  blocks  of  data  ( Cx,Cy,Cz,Ci,Cm,C„ )  for  6a  =  -18°  (ie,  configu¬ 
rations  3  -  4;  8  -  9),  and  three  incremental  blocks  for  S„  =  18°  (configurations 
2-  1;  5  —  4;  7  —  6).  The  result  of  this  differencing  shows  that  AC*,  A  Cy,  AC|  and 
A C„  are  close  to  zero  and  therefore  can  be  considered  independent  of  elevator  and 
rudder  control.  A Cz  and  A Cm  however,  do  show  evidence  of  coupling  but  when 
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compared  to  the  overall  magnitude  of  Cz  and  Cm,  the  Sa  increments  are  of  less  sig¬ 
nificance  than  the  influence  of  S,  and  £r.  This  is  convenient  because  now  the  aileron 
contribution  for  the  forces  and  moments  can  be  simply  added  as  an  increment  to 
the  elevator  and  rudder  data  blocks  as  follows: 

Cx  =  Cx{S„6rJ,a,  +  dCx{6aJ,a, 


Note  that  this  implies  that  when  Sa  =  0,  dCx  =  0,  etc. 

Since  only  one  incremental  data  block  is  needed  for  each  aileron  control  setting 
the  two  Sa  =  -18°  data  blocks  have  been  averaged.  Similarly  the  three  Sa  =  18° 
data  blocks  have  been  averaged.  Obviously  for  Sa  =  0°  all  incremental  data  is  zero 
and  no  data  block  is  needed.  This  new  data  base  scheme  is  shown  in  Figure  2. 

The  creation  of  the  incremental  aileron  data  blocks  also  solved  another  problem 
caused  by  limited  control  deflection  data.  The  following  control  deflections  were 
not  covered  in  the  NASA  LRC  test  sequence: 

1.  Sc  =  23°,  Sr  =  -25°,  Sa  =  0° 

2.  6,  =  -25°,  Sr  =  25°,  Sa  =  0° 

However,  the  following  control  deflections  with  aileron  were  available: 

3.  Se  =  23°,  Sr  =  -25°,  Sa  =  18° 

4.  S,  =  -25°,  6r  =  25°,  Sa  =  -18° 

By  subtracting  the  incremental  aileron  data  blocks  formed  above  from  the  18°  and 
-18°  aileron  configuration,  (3.  and  4.  above.),  data  for  zero  aileron  cases,  (1.  and 
2.  above.),  can  be  created. 

3.3.2  Generating  Positive  Sideslip  Angle  Data 

The  only  configuration  tested  at  NASA  LRC  for  both  positive  and  negative  d  angles 
was  for  St  =  Sr  =  Sa  =  0°,  where  the  /S  settings  were  -10°,0°,10°.  All  other  con¬ 
figurations  were  tested  for  f)  =  -10°  and  0°  since  preliminary  calculations  showed 
that  equilibrium  sideslip  angle  for  positive  spin  rates  could  be  expected  to  lie  within 
this  range.  For  the  general  unsteady  manoeuvres  occurring  during  spin  entry  and 
recovery  both  positive  and  negative  sideslip  angles  can  be  expected  to  occur.  How¬ 
ever  the  jd  =  10°  data  can  be  generated  by  assuming  that  the  aerodynamic  changes 
due  to  /9,  fl,  and  Sr  are  symmetrical,  although  in  practice  it  is  difficult  to  achieve 
symmetrical  conditions  in  a  wind-tunnel  environment.  Take  for  example  an  aircraft 
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spinning  with  the  following  conditions: 


8. 

- 

25° 

Sr 

= 

-25° 

8a 

= 

0° 

0 

= 

-10° 

a 

— 

60° 

Ob 

W 

= 

-0.6 

To  make  the  spin  rate  opposite  to  the  current  spin  rate,  (ie.  port  spin  to  starboard 
spin),  the  conditions  would  have  to  be  as  follows: 

6,  =  25° 

Sr  =  25° 

8a  =  0° 

0  =  10° 
a  =  60° 

$  =  0.6 

Note  that  8r,0,  and  only  have  changed  sign.  This  means  that  data  points  at 
negative  0  can  be  used  for  positive  0  with  opposite  spin  through  the  assumption 
of  spin  symmetry.  Below  is  an  example  of  how  the  positive  0  data  is  calculated. 
Figures  3  and  4  show  this  graphically. 


Cj{6e  =  25°, 5r  =  — 25° ,  /9  =  10°}  =  -Ct{8e  =  25°, «r  =  25°, /3  =  -10°}  (3) 

$  =  -0.6  —  0.6  $  =  0.6  -0.6 

a  =  0°  — *  90°  a  =  0°  -*  90° 

and  similarly  for  Cy,C„. 

Cm{£,  =  25°, 5r  =  -25°, 0  =  10q}  =  Cm{£,  =  25°, 6r  =  Zo°,0  =  -10°}  (4) 

$  =  -0.6  —  0.6  $  =  0.6  -0.6 

a  =  0°  -»  90°  a  =  0°  90° 

and  similarly  for  Cx,Cz- 

The  same  procedure  could  also  be  applied  to  the  aileron  incremental  block  data, 
but  this  was  not  required.  The  aileron  data  only  shows  minimal  change  due  to  0 
and  has  therefore  been  taken  to  be  independent  of  0.  All  aileron  incremental  data 
blocks  are  those  for  0  =  0°  and  are  constructed  as  follows: 
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dCx  =  fn{6.,a,ty)  (5) 

and  similarly  for  dCy  ,dCz}  dCi,  dCmi  dCn.  The  total  collated  data  for  the  model 
then  becomes: 


Cx 

=  Cx{6e,8,,0,a, 

nb\ 
2  Vi 

Cy 

=  Cy{8„8r,P,ct, 

Clb  \ 

W  j 

Cz 

=  Cz{8,,8r,0,a, 

Clb\ 

W! 

C, 

=  Ci{8„8,,f),a, 

nb\ 
2  V  s 

Cm 

=  8r,@,  a, 

06  \ 
2  VJ 

Cn 

=  Gn{8e,6r,fi,at, 

n&  \ 

w) 

+  dCx{8a,a,§} 

+  dCy{8a,a,$} 

+  dCz{Sa,a,  ^7} 

+  dCi{6a,a,  (6) 

+  dCm{8a,a,§} 

+  dCn{6a,a, 


This  is  depicted  in  Figure  5.  The  next  stage  iB  to  determine  a  means  by  which  to 
access  this  data  efficiently. 

3.3.3  Data  Base  Processing 

A  data  base  processing  program  library  has  been  developed  to  read,  store,  and  inter¬ 
polate/extrapolate  the  Wamira  spin  data  as  represented  in  the  data  base  schematic 
of  Figure  5.  Basically  each  force  and  moment,  (or  arbitrary  data),  has  been  defined 
as  a  variable  function  with  a  set  of  input  variables.  In  FORTRAN  the  axial  force 
coefficient  would  be  accessed  as  follows: 

CX  =  CX.DAT  A{DELE,  DELR,  BETA,  ALP  HA, OMEGA)  + 
DCX(DELA,  ALPHA, OMEGA) 

where  DELE,  DELR,  DELA,  BETA,  ALPHA,  OMEGA  can  be  any  arbitrary  value. 
All  interpolation/extrapolation  is  calculated  within  the  function  CX-DATA  and 
therefore  does  not  interfere  with  the  flow  of  the  equations  of  motion  for  the  aero¬ 
dynamic  routines. 

4  TRIM  ROUTINES 

It  is  preferable  to  begin  simulation  with  all  the  equations  of  motion  in  trim.  This 
means  that  all  accelerations  must  be  zero  along  and  about  the  body  axes  of  the 
aircraft.  The  Wamira  flight  dynamic  model  has  been  developed  to  trim  in  either 
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steady  spin  or  steady  rectilinear  flight  through  the  use  of  two  trimming  routines. 
Each  of  the  two  trim  conditions  is  calculated  in  two  trim  phases;  the  first  phase 
provides  an  approximate  initial  trim,  and  the  second  uses  these  values  in  a  more 
intensive  search  routine  to  provide  more  accurate  trim  conditions  using  the  flight 
dynamic  equations  within  the  ACSL  program. 

4.1  Steady  Spin  Trimming 

The  initial  approximate  spin  trimming  is  achieved  using  approximate  equilibrium 
spin  equations  and  follows  the  approach  presented  in  Reference  7.  Following  the 
assumptions  given  in  Reference  7,  Figure  6  shows  the  orientation  of  the  forces  and 
axes  of  an  aircraft  in  a  steady  spin.  The  equations  given  below  show  the  balance  of 
forces  in  a  steady  spin: 


mg  =  ±pV*SCD 

(7) 

mRsns  =  ~pViSCL 

(8) 

Note  that  in  this  approximation  the  side  forces  are  assumed  to  be  zero.  The  aero¬ 
dynamic  moments  of  Equation  (6)  are  balanced  by  the  following  inertial  moments: 

Cm 

=  (4/pSd6J)(!$),(/v  -  I z)  sin  2a cos* [a  +  0) 

(9) 

Ci 

=  (4/pS63)(^)*(/z  -  ly )  sin  at  sin  2{a  +  0) 

(10) 

Cn 

=  (4/pS6s)(^),(Jr  -  lx)  cos  a  Bin  2(<r  +  0) 

(11) 

where 


a 


These  equations  must  be  solved  simultaneously  to  trim  for  a  steady  spin.  A  nu¬ 
merical  solution  which  exploits  the  nature  of  the  forces  and  gives  added  physical 
insight  can  be  achieved  by  sequentially  processing  the  pitching,  rolling,  and  yawing 
moment  equations  (Equations  (9), (10), (11)  respectively).  For  constant  control  de¬ 
flections  the  inertia  terms  on  the  right  hand  side  of  the  equations  are  a  function  of 
a,  $,  and  0.  This  is  also  true  of  the  aerodynamic  coefficients  as  indicated  by  Equa¬ 
tion  (6)  of  Section  3.3.2.  The  a,  $  and  0  equilibrium  values  can  be  determined  as 
follows. 


tan 


-1 


V 

-1  r  pSbCL  , 

tan 


(12) 
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The  pitching  moment  equation  is  first  solved  for  $  for  the  range  of  test  a  values 
(here  0  — *  90°)  using  Equation  (9),  with  0  taken  as  zero  since  Cm  is  not  a  strong 
function  of  0  for  most  aircraft,  aerodynamically  and  inertially.  This  produces  a 
series  of  a  and  $  values  that  satisfy  the  pitching  moment  equation  as  shown 
in  Figure  7.  These  are  then  used  in  the  rolling  moment  equation  to  solve  for  0 
(Equation  (10)).  The  rolling  moment  equation  is  a  strong  function  of  0  in  both 
aerodynamic  and  inertial  terms, and  by  solving  for  0  with  the  previously  formed 
series  of  a  and  a  new  series  of  data  is  formed  that  satisfies  both  the  pitching 
and  rolling  moment  equations  as  shown  in  Figure  8. 

At  this  stage  if  required,  a  feed  back  loop  can  be  implemented  to  update  the  0  = 
0  condition  used  in  the  pitching  moment  equation  (Equation  (9))  to  the  new  values 
just  formed  in  Figure  8.  However,  this  is  only  required  if  the  aircraft  configuration 
being  studied  shows  signs  of  pitch  sensitivity  due  to  0  using  Equation  (9).  For  the 
Wamira  spin  model  data  the  correction  is  not  significant  and  does  not  require  the 
extra  calculation. 

Finally  the  yawing  moment  equation  (Equation  (11))  is  used  to  find  the  com¬ 
bination  of  a,  0  of  Figure  8  that  satisfy  all  three  moment  equation.  Each 
combination  is  substituted  into  both  the  aerodynamic  and  inertia  terms  of  the  yaw¬ 
ing  moment  equation  to  find  the  intersection  points  as  illustrated  in  Figure  9. 

The  intersection  of  the  inertial  and  aerodynamic  curves  identifies  the  combi¬ 
nation/s  of  a,  0  that  satisfy  all  three  moment  equations  for  equilibrium  in 
the  spin.  Only  those  intersections  where  the  aerodynamic  contribution  has  negar 
tive  slope  with  respect  to  the  inertia  curve  can  result  in  a  stable  spin,  (ie.  stable 
equilibrium).  Once  o,  and  0  have  been  found  for  equilibrium,  the  drag  and 
lift  equations  (Equations  (7)  and  (8)  respectively)  can  be  used  to  find  the  descent 
velocity  V  and  the  spin  radius  R$. 

The  initial  body  axis  rotation  rates  can  be  calculated  from  Reference  7  using 
the  helix  angle  from  Equation  (12)  as  follows: 

p  =  O  cos  o  cob(0  +  o}  ( 13) 

q  =  n8in(/S  +  a)  (14) 

r  =  fi  sin  a  cos(£  +  a)  (15) 

and  with  the  flight  path  angle  7  approximately  equal  to  90°,  9  and  <j>  are  calculated 
as  follows: 


9  ra  a  -  7  {7  <s  90°}  (16) 

<t>  ta  0  +  a  (17) 

This  trim  procedure  is  only  the  first  part  of  the  trim  sequence.  The  flight 
dynamic  spin  model  program  uses  these  initial  trim  results  to  perform  its  own  trim¬ 
ming  of  the  six  degree  of  freedom  flight  dynamic  equations  using  the  state  variables 
defined  in  Section  2.  This  trim  routine  is  called  POWIT,  (Reference  8),  and  is 


13 


a  multivariable  search  technique  for  the  solution  of  simultaneous  non-linear  equa¬ 
tions.  Basically  POWIT,  satisfies  a  given  a  set  of  condition  constraints,  by  solving 
the  full  aerodynamic  spin  equations  using  a  set  of  specified  trimming  variables  that 
are  incrementally  changed.  Solutions  are  found  using  a  combination  of  search  and 
gradient  techniques,  (detailed  descriptions  of  these  methods  can  be  found  in  Refer¬ 
ence  9).  Upon  exiting  this  routine  the  trimming  variables  define  an  accurate  trim 
state.  Figure  10  illustrates  the  procedure  with  a  block  diagram.  For  steady  spin, 
the  trimming  variables  and  condition  constraints  are  as  follows: 


'Prim  variables 
a 

P 

V 

e 

<t> 

p 

q 


Condition  constraints 
a  =  0 

4  =  0 

V  =  0 
9  =  0 

4  =  0 

P  =  ° 

?  =  0 
r  =  0 


tp  =  0  fixed  4  =  O 

for  given  control  setting  angles. 


4.2  Steady  Rectilinear  Flight  Trimming 

For  trimming  in  rectilinear  flight  the  initial  trim  estimate  is  supplied  by  manual  in¬ 
put.  The  initial  estimate  does  not  need  to  be  accurate  for  rectilinear  flight  because 
with  fewer  variables,  the  trimming  is  more  stable  than  for  the  spin.  Steady  rectilin¬ 
ear  trimming  is  performed  solely  using  POWIT  and  the  following  trimming  variables 
and  condition  constraints  are  required  to  satisfy  the  flight  dynamic  equation  block 
shown  in  Figure  10: 

Trim  variables 
a 
9 
4> 

6. 

Sr 

s* 

for  given  V  and  Thrust 

Note,  unlike  the  arrangement  for  trimming  in  the  steady  spin,  the  control  deflec¬ 
tions  are  not  fixed  but  are  used  instead  to  trim  the  aircraft  for  a  constant  velocity 


Conditions  constraints 
d  =  0 
V  =0 

4  =  o 
j  =  0 
r  =  0 

p  =  0 
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V  and  powerplant  thrust.  Since  the  rotary  balance  measurements  contain  asymme¬ 
tries  in  the  lateral  and  directional  forces  for  the  steady  symmetrical  Sight  conditions, 
small  control  deSections  Sa,  Sr  are  required  for  trim. 


5  DECOMPOSITION  OF  AIRCRAFT  ANGULAR 
MOTION 

The  basis  of  the  proposed  model  formulation  is  the  separation  of  the  total  angular 
motion  into  a  component  aligned  with  the  Sight  path  vector  and  a  residual  dis¬ 
turbance  component.  The  aerodynamic  forces  and  moments  due  to  the  former  are 
identically  those  measured  on  a  rotary  balance,  while  the  moments  due  to  the  resid¬ 
ual  disturbance  are  to  be  calculated  from  a  simple  aerodynamic  model.  The  net 
forces  at  the  c.g.  due  to  the  residual  component  have  been  taken  to  be  zero  in  the 
present  formulation.  In  other  words,  the  residual  component  is  modelled  simply  as 
a  redistribution  of  the  aerodynamic  load  without  any  change  in  the  resultant  X,Y, 
and  Z  components. 

The  method  of  separating  the  aligned  component  and  the  residual  disturbance 
component  was  introduced  by  Reference  10.  However,  the  equations  derived  assume 
the  aircraft  is  descending  vertically.  The  equations  presented  below  do  not  make 
this  assumption  and  therefore  allow  the  prediction  of  aircraft  behaviour  during  the 
transition  from  rectilinear  flight  to  steady  spin. 

The  rotation  rate  about  the  flight  path  vector  can  be  calculated  as  follows  (Ref¬ 
erence  11  Equation  (5.8,5)(c)). 

pw  =  pcoeacoBf)  (q  -  a)  sin  +  r  sin  a  cos  ^  (18) 

Note,  pw  is  equal  to  the  steady  rotation  rate  (fi)  which  is  used  for  the  rotary  balance 
test.  The  a  term  in  Equation  (18)  presents  a  problem  because  it  is  calculated  using 
the  forces  read  from  the  RB  data  base,  but  the  forces  can  only  be  calculated  after 
pw  is  found.  One  solution  is  to  iterate  the  system  of  equations  until  pw  and  a 
converge  on  a  steady  value.  Five  iterations  were  found  to  be  enough  for  an  accurate 
convergence  and  once  pw  (H)  has  been  calculated,  it  is  used  to  calculate  the  rotary 
and  disturbance  terms  of  p,  q,  r  which  are  derived  as  follows. 

From  Equation  (5.2,12)  of  Reference  11  p,  q,  and  r  may  be  calculated  using 
pw,  qw,  and  rw  through  rotations  a  and 


P 

Pw 

'  0  ' 

'  0  ‘ 

Lwb 

? 

- 

gw 

=  - 

0 

P  +  Lwb 

1 

r 

.  ** 

1 

0 
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(19) 


where  Lwb  is  the  matrix  for  transforming  vectors  from  body  axes  to  wind  axes  as 
follows: 


Lwb 


cos  a  cos  9  sin  /9 
—  coeasin^  cos ,9 
-  sin  a  0 


sin  a  cos  9 

-  sin  oi  sin  9 
cos  a 


(20) 


Now  since 


then 


Lbw  *  IwB 


(21) 


p 

q 

-  LBw 

pw 

qw 

r 

rw 

'  0  ' 

'  o  ' 

Lbw 

0 

1 

9  + 

l 

0 

(22) 


Equation  (4.5,5)  of  Reference  11  gives: 


Lbw 


cos  oi  cos  9 
sin  9 

sin  a  cos  9 


-cosasin9  -  sin ct 
cos  9  0 

-  sin  a  sin  9  cosat 


(23) 


By  substituting  Equation  (23)  into  Equation  (22)  p,  q,  and  r  can  be  separated  into 
components  due  to  the  rotation  about  the  wind  axes  (ie.  into  its  rotary  part)  and 
into  the  disturbance  components  as  follows: 


P 

1 

T 


Pw  cos  a  cos  9  |  -qw  cos  a  sin  9  ~  nv  sin  a  +  f)ain  a 
pw  sin  9  I  qw  cos  9  +  a 

pw  sinacos9  |  -qw  sin  a  sin  9  +  rw  cos  a  -  $  cos  a 


rotary 


disturbance 


(24) 


Now  since  the  state  equations  of  the  Sight  dynamic  model  calculate  p,  q,  and  r 
directly,  the  disturbance  terms  can  be  evaluated  as  follows  to  reduce  the  number  of 
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mathematical  calculations 


Po  =  P~Pr 

1o  =  q  ~  qr 

r0  =  r  -  rr 


(25) 


Once  the  disturbance  rotation  rates  are  calculated  using  Equations  (24)  and  (25) 
the  disturbance  moments  can  be  found  using  the  following  equation  and  the  calcu¬ 
lated  data  base: 


Ci. 

Cip  cu  Ci, 

Po 

Cm0 

= 

Cmp  Cm^  Cra, 

7o 

Cn9 

.  Cn„  C„r 

Only  the  diagonal  terms  of  the  moment  coefficient  matrix  have  so  far  been  used  to 
calculate  the  disturbance  moments  as  these  terms  appear  to  be  the  most  significant. 
The  effect  of  the  other  terms  has  been  assumed  to  be  small  and  they  are  taken  as 
zero.  Future  investigations  into  the  disturbance  coefficients  will  be  undertaken  to 
evaluate  the  influence  of  these  smaller  terms  on  predicted  spin  behaviour. 

5.1  Small  Disturbance  Data 

The  rotary  balance  data  was  measured  with  the  Wamira  spin  model  performing 
steady  rotations  about  an  axis  aligned  with  the  freestream.  This  provides  the 
steady  aerodynamic  moment  data  but  does  not  supply  the  disturbance  data  that  is 
required  in  this  model  formulation  to  represent  the  forces  occurring  in  non-steady 
flight.  The  small  disturbance  coefficients  included  in  Equation  26  represent  the 
aerodynamic  moments  due  to  body  rotation  rates.  The  main  diagonal  terms,  which 
have  been  included  in  the  model,  generally  contribute  damping  forces  in  the  aircraft 
response.  A  simple  line  vortex  model  has  been  developed  to  calculate  these  terms. 

The  line  vortex  model  estimates  the  aerodynamic  coefficients  during  steady  ro¬ 
tation  using  line  vortex  representations  of  wing,  tailplane,  and  fin  lift.  The  vortex 
model  is  matched  to  the  measured  rotary  balance  data,  and  once  a  satisfactory 
match  has  been  estimated,  small  disturbances  in  angular  body  rates  are  applied 
and  associated  changes  in  aerodynamic  coefficients  are  calculated  to  provide  distur¬ 
bance  moment  derivatives  in  Equation  26.  The  disturbance  moment  coefficients  are 
calculated  over  the  full  range  of  a  and  ^  to  create  the  necessary  data  base.  The 
effect  of  0  is  neglected.  Before  this  data  base  can  be  accessed  the  aircraft  rotary 
and  disturbance  rotation  rates  must  be  calculated.  A  report  on  this  computational 
procedure  is  in  preparation. 
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6  EXAMPLE  APPLICATION  OF  THE  SPIN  MODEL 

6.1  Equilibrium  Spin  Prediction 

The  flight  dynamic  spin  model  can  be  used  to  calculate  possible  spin  modes  through 
the  use  of  the  trimming  routines  explained  in  Section  4.1.  Table  1  shows  a  compar¬ 
ison  of  a  spin  mode  given  in  Reference  5,  calculated  by  the  same  technique  used  by 
the  present  model  for  the  initial  trim,  and  the  same  spin  mode  calculated  with  the 
flight  dynamic  spin  model.  Other  modes  can  be  identified  by  changing  the  setting  of 
the  control  deflections.  Note  the  difference  between  the  spin  prediction  of  FDSM(l) 
and  FDSM(2)  of  Table  1.  The  spin  trim  equations  of  Section  4.1  are  represented  by 
FDSM(l)  and  assume  that  all  side  forces  are  zero.  FDSM(2)  however,  represents 
the  full  flight  dynamic  spin  equations  of  motion  and  does  include  side  force  effects. 


Table  1:  Spin  Prediction  Comparison 


a 

■!!■ 

WmvMtk 

mm 

0 

mm 

Ref.  5 

o 

0.28 

o 

-4 

MEM 

FDSM(1Y 

mn 

0.276 

gflTil 

-3.58 

EEH 

wmsmera i 

gEEl 

0.341 

mm 

-5.14 

KE3 

(for  5«  =  -25°,  6,  =  -25°,  Sa  =  0°) 


*  Initial  »pin  approximation 
**  POWJT  6  d.o.r  trim 


6.2  Spin  Manoeuvre 

To  show  the  full  capability  of  the  flight  dynamic  spin  model  a  manoeuvre  that  covers 
both  rectilinear  and  spinning  flight  will  be  presented  as  an  example.  Figures  11-16 
show  the  time  history  traces  for  the  manoeuvre  which  is  constructed  as  follows. 

Firstly,  the  flight  dynamic  spin  model  is  trimmed  in  rectilinear  flight  for  a  flight 
path  velocity  of  90  kts  and  at  an  altitude  of  10000  ft.  The  trim  results  are  a  =  8.2°, 
•y  =  -6.7°  with  the  powerplant  thrust  set  to  be  zero,  and  the  control  trim  values 
6,  =  -2.6°,  6r  =  1.8°,  and  Sa  =  -1.9°.  Note  the  offset  in  fy  and  6„  required  to 
balance  non-zero  lateral  and  directional  moments  measured  in  the  RB  test  to  keep 
the  aircraft  in  a  wing  level  descent. 

With  the  aircraft  in  steady  descent  for  5  seconds  the  elevator  is  moved  t of,  = 
-25°  at  time  t  =  5  seconds.  The  aircraft  pitches  up  and  stalls  at  about  t  =  6 
seconds  at  which  time  full  starboard  rudder  is  applied  (Sr  =  -25°)  to  yaw  the 
aircraft  quickly  and  promote  a  fast  spin  entry.  The  controls  remain  fixed  until  time 
t  —  45  seconds  when  full  anti-spin  rudder  is  used  (S,  —  25°)  to  recover  from  the 
spin.  After  another  3  seconds  the  aircraft  has  recovered  but  still  slowly  rolls  until 
at  time  t  =  48.5  seconds  the  aircraft  begins  to  roll  in  the  opposite  direction.  At  this 
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instant  the  elevator  and  rudder  controls  are  brought  to  the  initial  trim  values.  The 
aircraft  then  dives,  buiids  up  speed  and  lift,  then  rounds  out  to  almost  level  Sight 
with  a  heading  45°  starboard  from  where  it  started  and  has  lost  almost  9000  ft 
altitude.  All  these  control  deflections  are  shown  in  Figure  11. 

Figure  12  shows  the  aircraft  decreasing  from  90  to  75  kts  during  the  spin  entry, 
while  a  converges  to  48“  (expected  from  Table  1)  and  -y  approaches  -90°  (vertical 
descent).  The  velocity  trace  also  shows  the  sharp  increase  during  the  recovery 
(almost  225  lets)  and  then  a  decrease  as  lift  increases  in  the  pull-out.  Note  that 
during  the  recovery  stage  (time  t  >  48.5  seconds)  that  the  high  spin  a  has  returned 
to  standard  flight  values  of  about  3°  to  5°, 

The  time  history  traces  of  4,  9,  t/i  are  shown  in  Figure  13.  During  the  Bpin 
the  spin  rate  converges  to  ^  =  0.36,  then  returns  to  zero  at  recovery.  The  oscil¬ 
latory  motion  displayed  during  spin  entry  and  recovery  is  due  to  a  lowly  damped 
Dutch  roll  mode.  A  more  detailed  investigation  into  the  aerodynamic  data  at  low 
a  will  be  reported  later.  Figure  14  shows  the  body  rotation  rates  and  normal  accel¬ 
eration  during  the  manoeuvre.  The  pitch  rate  trace  illustrates  the  sharp  pitching 
due  to  the  elevator  control  deflections  at  time  t  —  5  seconds  for  the  spin  entry  and 
at  time  t  =  48.5  seconds  during  the  recovery.  The  yaw  rate  trace  displays  the  high 
continuous  yaw  rates  that  can  only  be  achieved  during  a  spin  (about  2.1  rad/sec). 
The  dashed  lines  in  the  angular  rotation  rate  manoeuvre  traces  represent  the  angu¬ 
lar  velocity  components  relative  to  axes  aligned  with  the  flight  path  vector  while  the 
solid  lines  represent  the  total  angular  rates.  The  difference  between  these  two  curves 
therefore  represents  the  disturbance  component  of  angular  rotation  rate  given  by 
the  equations  in  Section  5.  The  difference  between  these  curves  for  steady  pitch 
rate  arises  because  the  RB  contribution  is  related  to  the  sideslip  at  zero  spin  radius 
(Equation  (24))  while  the  total  pitch  rate  is  related  to  sideslip  and  to  the  helix  angle 
associated  with  the  spin-radius  (see  approximate  Equation  (14)).  Aircraft  position 
is  shown  in  Figure  15  against  time  while  Figure  16  plots  altitude  and  displacement 
east  (Y)  against  displacement  north  (X). 

The  time  history  traces  provide  a  large  amount  of  detailed  information  but  do 
not  quickly  illustrate  the  flight  path  and  aircraft  attitudes  during  the  manoeuvre. 
A  general  application  graphics  routine  has  been  developed  to  read  data,  (not  nec¬ 
essarily  time  history),  and  plot  it  in  three  dimensions  (3D).  Basically  the  graphics 
program  reads  a  time  history  then  plots  the  flight  path  on  to  which  an  aircraft 
image  is  overlaid  at  constant  time  intervals  with  full  directional  orientation. 

The  previous  manoeuvre  has  been  graphically  illustrated  in  this  way  and  is 
shown  in  Figures  17  -  20.  The  full  manoeuvre  is  plotted  in  Figure  17  with  an 
aircraft  displayed  every  2.5  seconds.  The  time  interval  must  differ  from  the  spin 
period  which  in  this  case  is  2.3  sec/turn,  or  else  the  same  view  of  the  aircraft  would 
be  seen  repeatedly. 

The  3D  plot  now  makes  it  easier  to  visualise  the  attitudes  of  the  aircraft  from 
trim  to  stall,  spin  entry  to  spin,  and  from  spin  to  recovery.  Figure  18  displays  a 
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more  detailed  view  of  the  spin  entry  and  shows  how  the  aircraft  pitches  up  into  a 

9taU  then  banks  starboard  as  the  rudder  is  applied.  Figure  19  shows  the  <ur«aft 
staU  then  banks  to  star  detai]8  the  recovery  stage.  The  black 

7  CONCLUSION 

‘.zxs£s~ rcSHSarS 

datafromalmevortexmodd.  o(  the  progrMn  to  simulate  spin 

be  reported  later.  The  mam  aspects  are: 

a)  Details  of  the  line  vortex  model  used  to  calculate  the  small  disturbance  data, 

b)  Effects  of  the  cros^terms  in  the  small  disturbance  coefficient  matrix, 

c)  The  specific  cause  of  the  predicted  oscillatory  behaviour  at  about  spin  fre¬ 
quency  during  spin  entry  and  recovery, 

d)  Validation  of  the  modelling  by  reference  to  dynamic  testinj >  “  * y^aut 
model  conducted  at  the  NASA  LRC.  Full  k*  data  are  not 

the  Wamira  project  did  not  proceed  to  the  flight  prototype  stage.  Refer 
discusses  the  model  application  to  spin  recovery. 
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Figure  10.  Flight  dynamic  spin  model  block  diagram. 
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Figure  14.  Manoeuvre  time  history  traces, 
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